AD-A147  406  CALCULATIONS  OF  HIGH-CURRENT  CHARACTERISTICS  OF  SILICON  i/i 
DIODES  AT  MICROWAVE  FREQUENCIES^)  HARRV  DIAMOND  LABS 
ADELPHI  MD  A  L  WARD  OCT  84  HDL-TR-2057 


UNCLASSIFIED 


F/G  9/1 


NL 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  rWbon  Dots  F.nfrod) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

HDL-TR-2057 


4.  TITLE  (and  Subtitle) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3  RECIPIENT'S  CATALOG  NUMBER 


5  TYPE  OF  REPORT  b  PERIOD  COVERED 


Calculations  of  High-Current  Characteristics  of 
Silicon  Diodes  at  Microwave  Frequencies 


7.  AUTHOR!*,) 

Alford  L.  Ward 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U  S.  Army  Materiel  Command 
5001  Eisenhower  Ave 

Alexandria.  VA  22333-0001  _ 


4.  MONITORING  AGENCY  NAME  a  ADDRESS (If  dlffarant  I 


»6.  DISTRIBUTION  STATEMENT  (of  thla  Raport) 


Technical  Report 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBER( a) 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  6  WORK  UNIT  NUMBERS 

Program  Ele:  621 20A 
DA  Project:  1L162120AH25 


12.  REPORT  OATE 

October  1984 


13.  NUMBER  OF  PAGES 

25 


i  Controlling  Otflca)  15.  SECURITY  CLASS,  (of  thla  raport ) 

UNCLASSIFIED 


ts«.  OECLASSI  FI  CATION/ DOWN  GRADING 
SCHEOULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  fh#  abatract  antarad  In  Block  20,  II  dlffarant  from  Raport) 


18.  SUPPLEMENTARY  NOTES 


HDL  Project:  E013E1 
AMS  Code:  612120H25001 


19.  KEY  WORDS  (Contlnua  on  ravaraa  alda  It  nacaaaary  and  Idantlfy  by  block  numbar) 

Silicon  diode 
Microwave  frequency 
Second  breakdown 
Forward  conduction 
Avalanche  multiplication 


ZO.  ABSTRACT  fCoottmio  as  r«WM  at*  ft  nmconmr  Idoiultr  ftr  block  numbot) 

Dynamic  current/voltage  characteristics  are  calculated  for  silicon  diodes  at  microwave  frequencies. 
Forward  and  reverse  half  cycles  are  studied  separately  Little  current  flows  unless  the  forward  half  period 
exceeds  the  carrier  transit  time:  then,  as  the  frequency  decreases,  the  peak  forward  current  increases 
essentially  linearly  with  the  half  period,  until  the  dc  value  is  attained.  Avalanche  multiplication  is  noted  in 
the  forward  direction  at  about  one  half  of  the  reverse  breakdown  voltage  Below  the  reverse  breakdown 
voltage,  the  current  consists  solely  of  the  stored  forward  charge  Above  the  reverse  breakdown  voltage, 
the  peak  current  has  a  maximum  at  a  frequency  which  depends  on  the  diode  width  and  temperature.  > 


EDITION  OF  '  NOV  SS  IS  OBSOLETE 


UNCLASSIFIED 

security  classification  of  THIS  PAGE 


i  Data  Bn  tar  ad) 


r 


r 


u 


CONTENTS 


1 .  INTRODUCTION 


2.  COMPUTER  PROGRAM 


3.  FORWARD  HALF  CYCLE 


4.  REVERSE  HALF  CYCLE 


5.  REVERSE  AVALANCHING  .  10 

6.  FORWARD  AVALANCHING  .  13 

7.  TEMPERATURE  DEPENDENCE  .  14 

8.  CONCLUSION  .  18 

LITERATURE  CITED  .  19 

DISTRIBUTION  .  21 


FIGURES 


1.  External  circuit  used  in  calculations 


2.  Forward  cur rent /voltage  characteristics  for  1-pm  p-i-n  diode  for 
constant  voltage  and  various  frequencies  . 


3.  Peak  forward  current  density  as  a  function  of  applied  frequency 
for  various  diode  widths  and  voltages  . 


4.  Current,  voltage,  and  power  during  forward  and  reverse  half  cycles 
with  an  applied  voltage  of  60  V  peak  on  a  4-ym  diode  . 


5.  Maximum  current  density  as  a  function  of  frequency  for  reverse-bias 
(half-cycle)  0. 2-ym  diode  . . . 


6.  Dynamic  current/voltage  reverse  characteristics  for  100-V  peak 
applied  voltage  . . . . 


7.  Overvoltage  as  a  function  of  rate  of  voltage  rise  for  data  for 

figure  5  . . . 

8.  Comparison  of  forward  and  reverse  current  peaks  as  a  function  of 

frequency  for  0.2-ym  diode  . . 


9.  Forward  turn-on  transients  for  100-V  maximum  applied  voltage  at 

160  GHz  .  13 


VVV 


FIGURES  (Cont'd) 


1 


INTRODUCTION 


Second  breakdown  due  to  a  sinusoidal  voltage  waveform  is  of  interest  to 
the  electrical-overstress/electostatic-discharge  community  as  well  as  to  the 
microwave  and  electromagnetic  radiation  communities.  By  far  the  majority  of 
second -breakdown  studies  have  been  for  single -polarity  square  pulses.  Square- 
pulse  studies  are  convenient  for  both  the  experimentalist  and  the  theoreti¬ 
cian.  Several  investigators  have  predicted  second  breakdown  due  to  sinusoidal 
waveforms  using  suitable  modification  of  square-pulse  data. 

Whalen  and  Domingos1  have  compared  square-  and  rf -pulse  overstress  data 
for  bipolar  uhf  transistors.  They  reported  that  the  pulse  energy  required  to 
cause  failure  is  lowest  for  reverse-polarity  square  pulses,  intermediate  for 
rf  pulses  (240  MHz),  and  highest  for  forward-polarity  square  pulses.  Their 
finding  that  intermediate  energy  causes  failure  for  rf  pulses  is  logical  in 
the  thermal  model  of  second  breakdown,  even  though  the  thermal  model  itself 
does  not  distinguish  between  forward  and  reverse  breakdown.  Again,  the  ther¬ 
mal  model  is  inadequate  to  study  the  frequency  dependence  of  rf-pulse  second 
breakdown.  In  general,  the  thermal  model  is  valid  for  breakdown  times  greater 
than  the  thermal  time  constant  of  the  device,  usually  on  the  order  of  a  micro¬ 
second.  For  shorter  breakdown  times,  the  electrical  characteristics  dominate; 
this  regime  has  usually  been  called  the  current-mode2  second  breakdown. 
Greater  frequency  dependence  of  rf-pulse  second  breakdown  should  be  expected 
in  the  current-mode  time  regime. 

The  calculations  reported  herein  were  made  with  the  Harry  Diamond  Labora¬ 
tories  (HDL)  DIODE  computer  program  that  has  been  used  successfully  to  calcu¬ 
late  second  breakdown3  in  reverse-bias  junctions  and  to  study  various  modes  of 
avalanche  oscillations4  (IMPATT,  TRAPATT,  MULTIPATT*).  The  most  important 
limitations  of  the  present  program  are  (1)  it  is  limited  to  one-dimensional 
calculations,  (2)  it  neglects  thermal  conduction,  and  (3)  it  is  limited  to  one 
polarity  in  a  given  calculation. 

Rf  damage  of  a  semiconductor  junction  is  caused  by  the  temperature  rise  in 
the  high-field  junction  area.  The  disparity  of  the  electronic  time  scale 
(dielectric  time  constant  of  the  order  of  picoseconds)  and  the  thermal  time 


1J.  J.  Whalen  and  H.  Domingos,  Square  Pulse  and  RF  Pulse  Overstressing  of 
UHF  Transistors ,  Proceedings  of  the  Electrical  Overstress/Electrostatic 
Discharge  Symposium,  EOS-1  (1979),  140-146. 

2J.  H .  Fee,  W.  J.  Orvis,  L.  C.  Martin,  and  J.  C.  Peterson,  Modeling  of 
Current  and  Thermal  Mode  Second  Breakdown  Phenomena,  Proceedings  of  the 
Electrical  Overstress/Electrostatic  Discharge  Symposium,  EOS-4  (1982),  76-81. 

3 A.  L.  Ward,  Calculations  of  Second  Breakdown  in  Silicon,  Harry  Diamond 
Laboratories ,  HDL-TR-1978  (August  1982). 

4 A.  L.  Ward,  Modes  of  Avalanche  Oscillations  in  Silicon  Diodes,  IEEE  Trans. 
Electron  Devices,  ED- 2 5  (1978),  683-687. 

*  IMPATT— impact  avalanche  transit  time/  TRAPATT— trapped  plasma  avalanche 
transit  time;  MULTIPATT — multiple  avalanche  transit  time. 
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constant  (on  the  order  of  microseconds)  prevents  the  calculation  (at  a  reason¬ 
able  computer  cost)  of  a  sufficient  temperature  rise  to  cause  damage.  How¬ 
ever,  the  rate  of  temperature  rise  may  be  calculated  at  selected  temperatures; 
then  the  time  required  to  reach  the  damage  temperature  may  be  approximated. 

Damage  may  be  the  result  of  thermal  displacement  of  dopant  atoms  and/or  the 

melting  of  channels. 

This  investigation  studies  the  basic  dynamic  characteristics  of  silicon 
junction  diodes  upon  application  of  a  high-power  rf  pulse.  The  next  section 
describes  the  computer  program.  The  forward  half  cycle  and  the  reverse  half 
cycle  are  treated  separately  in  the  following  sections. 

2.  COMPUTER  PROGRAM 

A  computer  program  that  has  been  used  to  study  square-pulse  second  break¬ 
down  in  silicon  diodes3 *  has  now  been  used  to  study  the  microwave  characteris¬ 

tics  of  silicon  diodes. 

The  computer  program  was  designed  to  study  only  reverse-bias  breakdown, 
but  it  can  be  used  to  study  large  forward-bias  current  when  the  diffusion 
component  of  the  current  is  small  compared  to  the  conduction  component.5 
Negative  fields  tend  to  lead  to  calculation  instability;  therefore,  the  for¬ 
ward  and  reverse  half  cycles  of  the  microwave  signal  must  be  calculated  sepa¬ 
rately.  The  charge  distributions  at  the  end  of  the  forward  half  cycle  are 
used  as  initial  conditions  for  the  reverse  half  cycle.  Space -charge -free 
conditions  (i.e.,  when  the  mobile  charges  equal  the  fixed  charges)  are  used  as 
initial  conditions  for  the  forward  half  cycle. 

The  circuit  simulated  is  usually  chosen  to  consist  of  a  small  capacitance 
shunting  the  diode,  with  the  voltage  applied  through  a  series  resistance. 
This  circuit  is  shown  in  figure  1.  The  applied  voltage  may  be  chosen  to  be 
constant,  to  have  an  incremental  step  voltage,  or  to  have  the  applied  sinu¬ 
soidal  half  wave  used  mainly  in  this  study.  The  initial  voltage  need  not  be 
set  equal  to  the  applied  voltage.  Among  other  available  external  circuits  is 
one  with  an  inductance  in  series  with  the  external  resistance.  Calculations 
may  be  made  at  a  constant  temperature,  or  the  temperature  can  be  calculated  to 
increase  due  to  the  local  power  density.  Results  reported  in  sections  3 
through  6  are  for  room-temperature  calculations;  temperature  effects  are 
considered  in  section  7. 


3A.  L.  Ward,  Calculations  of  Second  Breakdown  in  Silicon,  Harry  Diamond 

Laboratories,  HDL-TR-1978  ("August  1982). 

SA.  L.  Ward,  Oscillating  Voltage  Pulses  and  Second  Breakdown,  Proceedings  of 
the  Electrical  Overstress/Electrostatic  Discharge  Symposium,  EOS-2  (1980), 
130-139. 
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Figure  1.  External  circuit  used  in 
calculations. 


3.  FORWARD  HALF  CYCLE 

We  first  discuss  calculations  of  the  characteristics  of  the  forward-bias 
half  cycle  of  diodes  subject  to  an  applied  microwave  voltage.  Figure  2  shows 
calculations  for  a  1-um  silicon  p-i-n  diode.  The  series  resistance  is  50  ft, 
the  shunt  capacitance  is  0.2  pF,  and  the  area  is  1  *  10“6  cm2.  The  low-level 
recombination  lifetime  is  50  us  in  the  intrinsic  material;  the  effective  life¬ 
time  decreases  as  the  current  increases.  Figure  2(a)  shows  the  current  densi¬ 
ty  as  a  function  of  time  for  four  applied  frequencies;  also  shown  is  the  turn¬ 
on  transient  for  a  fixed  voltage.  The  carrier  transit  time  across  a  1-um 
diode  at  saturated  carrier  velocities  is  10  ps.  Examples  of  charge  and  field 
distributions  during  a  turn-on  pulse  may  be  found  in  section  6.  Figure  2(a) 
shows  that  four  or  five  transit  times  are  required  for  the  peak  microwave 
current  pulse  to  approach  the  constant-voltage  saturation  current. 

Figure  2(b)  shows  the  diode  voltage  as  a  function  of  time  for  the  same 
applied  voltage.  Note  that  the  voltage  peaks  for  the  constant  applied  voltage 
case.  The  decrease  in  diode  voltage  is  due  to  the  current  flow  in  the  series 
resistance.  The  voltage  drop  is  greater  for  large  diode  areas  or  for  large 
series  resistances.  The  peaking  of  the  voltage  during  the  turn-on  transient 
has  been  observed  in  our  laboratories  and  by  others.  The  current/voltage 
characteristics  for  the  various  frequencies  are  shown  in  figure  2(c).  The 
counterclockwise  current/voltage  excursion  indicates  the  diode's  inductive 
reactance,  a  result  of  the  time  required  for  the  current  to  increase. 

The  current-density  maxima  are  plotted  as  a  function  of  frequency  in 
figure  3  for  the  1-um  p-i-n  diode  at  two  applied  voltages.  Also  shown  are 
similar  plots  for  4-,  20-,  and  100-pm  diodes. 

The  probable  usage  of  diodes  of  these  various  widths  can  be  determined  by 
the  reverse-breakdown  voltage.  The  breakdown  voltage,  VB,  as  a  function  of 
diode  width  for  silicon  diodes  which  are  just  punched  through  at  first  break¬ 
down,  is  given  by3 


where  the  diode  width,  d,  is  given  in  micrometers.  For  example,  1-um  diodes 
have  a  reverse  breakdown  of  about  21  V  while  100-pm  diodes  break  down  at  about 
1100  V.  The  cutoff  frequencies  shown  in  figure  3  are  the  upper  frequencies 
usable  for  each  diode  width. 


3 A.  L.  Ward,  Calculations  of  Second  Breakdown  in  Silicon,  Harry  Diamond 
Laboratories,  HDL-TR-1978  (August  1982 )• 
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Figure  2.  Forward  current/voltage  characteristics  for  1-um  p-i-n  diode  for 
constant  voltage  and  various  frequencies:  (a)  current  density  as  a  func¬ 
tion  of  time,  (b)  voltage  as  a  function  of  time,  and  (c)  current  density  as 
a  function  of  voltage. 

The  4-um  diode  of  figure  3  is  N-type,  doped  to  5  x  1015  cm-3,  whereas  the 
others  are  p-i-n  diodes.  Past  experience  has  shown  that  doping  level  has 
little  effect  on  the  diode  forward  characteristics  at  the  high  current  levels 
considered  here.  The  dc  current  levels  at  the  same  applied  voltage  are  indi¬ 
cated  at  the  left  margin  of  the  figure.  The  maximum  frequency  shown  for  each 
diode  is  approximately  the  inverse  of  the  transit  time  (at  saturated  velocity) 
for  each  device.  It  may  be  noted  that  more  transit  times  are  required  for 
current  saturation  for  the  longer  diodes.  With  this  in  mind,  we  find  that  our 
calculations  seem  compatible  with  the  measurements  of  Caulton  et  al.6  They 
measured  voltage  waveforms  for  diodes  of  about  100  to  300  pm  and  observed 
compression  of  the  voltage  waveform  due  to  forward  current  flow.  One  diode  of 
unspecified  thickness  showed  no  compression  above  41  MHz,  but  sharp  compres¬ 
sion  at  lower  frequencies.  They  found  that  this  frequency  increased  for 
decreased  diode  thickness,  in  agreement  with  our  calculations. 

6M.  Caulton,  A.  Rosen,  P.  J.  Stabile,  and  A.  Gombar,  p-i-n  Diodes  for  I«ow- 
Frequency  High-Power  Switching  Applications ,  IEEE  Trans.  Microwave  Theory 
Tech.,  MIT -30  (1982),  875-882. 
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Figure  3.  Peak  forward  current 
density  as  a  function  of  applied 
frequency  for  various  diode  widths 
and  voltages.  Diagonal  lines  are 
drawn  with  slope  of  -1 . 
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4.  REVERSE  HALF  CYCLE 


As  stated  earlier,  the  mobile  charge  distributions  at  the  end  of  the 
forward-bias  pulse  are  used  as  initial  conditions  for  the  reverse-bias 
pulse.  The  two  calculations  are  then  matched  at  the  zero  voltage  point  of  the 
applied  voltage  waveform.  This  entails  minor  extrapolations  of  the  calculated 
diode  voltage  and  current  toward  the  crossing  point.  The  diode  current  is 
zero  for  zero  diode  voltage  since  only  the  conduction  current  is  plotted. 
Diffusion  and  displacement  currents  are  calculated  but  are  not  plotted  here. 
Figure  4  shows  the  composite  waveform  for  a  4-ym  N-type  diode  doped  to  5  x 
1015  cm-3.  The  maximum  applied  voltage  is  60  V  for  each  bias  direction.  The 
series  resistance  is  again  50  ft,  the  shunt  capacitance  is  0.4  pF,  and  the  area 
is  1  x  10"5  cm2.  The  voltage  compression  in  the  forward  half  cycle  is  about 
1/3,  as  compared  to  the  reverse  half  cycle.  As  noted  earlier,  this  depends  on 
the  diode  area  used  as  well  as  the  series  resistance. 

The  calculated  power  dissipated  is  also  shown  in  figure  4.  Note  that 
little  power  is  dissipated  in  the  reverse  half  cycle  compared  to  the  forward 
half  cycle.  This  changes  drastically  when  the  reverse -breakdown  voltage  is 
exceeded.  The  shoulder  on  the  current  decay  (at  about  250  ps)  is  due  to  the 
transition  from  space-charge-limited  to  space-charge-free  conditions.5 

5A.  L.  Ward,  Oscillating  Voltage  Pulses  and  Second  Breakdown,  Proceedings  of 
the  Electrical  Overstress/Electrostatic  Discharge  Symposium,  EOS-2  (1980), 
130-139. 
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Figure  4.  Current,  voltage,  and 
power  during  forward  and  reverse 
half  cycles  with  an  applied  volt¬ 
age  of  60  V  peak  on  a  4-pm  diode. 


5.  REVERSE  AVALANCHING 


The  reverse  current  may  exceed  the  forward  current  when  the  avalanche 

voltage  is  exceeded.  Maximum  current  densities  calculated  as  a  function  of 

frequency  are  shown  in  figure  5  for  a  0.2-ym  p-i-n  diode.  This  diode  has  an 
avalanche  breakdown  voltage  of  11  V;  narrower  diodes  will  have  field-emission- 
assisted  breakdown  (not  calculable  with  the  present  program).  The  dropoff  of 
the  current  peak  at  high  frequencies  was  expected  (but  not  the  sharpness), 
although  that  at  lower  frequencies  was  not.  The  maximum  current  shown  for  an 

applied  voltage  of  100  V  at  the  lowest  frequency  is  lower  than  that  shown  for 

the  same  dc  voltage,  because  only  the  initial  peak  of  a  relaxation  oscil¬ 
lation4  is  shown.  Assuming  an  area  of  1  x  10”6  cm2,  appropriate  for  this 
diode,  a  maximum  instantaneous  power  of  over  1000  W  is  dissipated  for  the 
highest  applied  voltage.  This  heat  is  generated  in  a  volume  of  2  x  10”11  cm3, 
so  that  the  power  density  is  5  x  1013  W/cm3  and  the  temperature  rise  exceeds 
1  K  per  picosecond.  This  rate  of  rise  depends  only  on  the  current  density  and 
the  electric  field  and  is  thus  independent  of  the  area  chosen.  If  the  current 
were  constant  at  the  peak  current,  the  diode  would  burn  out  in  about  1  ns, 
making  an  experimental  measurement  of  pulse-length  dependence  very  diffi¬ 
cult.  The  dashed  curve  for  200  V  applied  is  calculated  with  a  shunt  capaci¬ 
tance  of  0.1  pF,  and  the  others  with  0.2  pF. 


The  dropoff  of  the  peak  current  as  the  frequency  decreases  may  be  under¬ 
stood  with  the  aid  of  figure  6,  where  the  current/voltage  dynamic  characteris¬ 
tics  are  shown  for  100  V  applied.  An  area  of  1  *  10-6  cm2  is  assumed.  The 
characteristics  are  traversed  in  a  counterclockwise  direction  with  time.  This 
is  a  result  of  the  inductive  effect  of  carrier  transport.  Figure  6  shows  that 
the  breakdown  (maximum)  voltage  increases  with  frequency,  and  the  current  peak 
increases  with  the  voltage  peak. 


L.  Ward,  Modes  of  Avalanche  Oscillations  in  Silicon  Diodes,  IEEE  Trans. 
Electron  Devices,  ED- 25  (1978),  683-687. 
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Figure  5.  Maximum  current  density  as 
a  function  of  frequency  for  reverse- 
bias  (half-cycle)  0.  2-ym  diode. 
Parameter  is  peak  applied  voltage. 


Figure  6.  Dynamic  current/voltage 
reverse  characteristics  for  1 00-V 
peak  applied  voltage.  Parameter 
is  angular  frequency,  w  =  2irf. 


Earlier  calculations3  have  shown  that  second-breakdown  overvoltages  in¬ 
crease  to  the  2/3  power  of  the  rate  of  rise  of  the  applied  voltage.  Figure  7 
is  a  plot  of  the  over-voltage  as  a  function  of  the  rate  of  voltage  rise  for 
the  characteristics  of  figure  6.  A  straight  line  results  on  this  log-log  plot 
when  a  breakdown  of  11  V  is  assumed;  deviations  result  from  assumed  breakdown 
voltages  of  10  or  12  V.  This  breakdown  voltage  agrees  with  that  found  in 
figure  6  when  the  current  peaks  are  extrapolated  to  zero  current.  The  static 
characteristic  indicates  a  differential  resistance  of  0.55  (2.  This  is  less 
than  the  forward  differential  resistance  in  the  same  current  range  <ind  may  be 
correlated  with  the  lower  power  for  second  breakdown  in  the  reverse-  as  com¬ 
pared  to  the  forward-bias  direction.  Further  calculations  will  be  required  to 
determine  whether  the  slope  of  0.58  seen  in  figure  6  applies  generally  to 
first-breakdown  voltage  as  compared  to  the  slope  of  0.67  for  second  breakdown. 

Comparisons  of  forward  and  reverse  current  peaks  as  a  function  of  frequen¬ 
cy  are  shown  in  figure  8  at  two  applied  voltages.  The  sharp  high-frequency 
cutoff  of  the  reverse-bias  current  peak  is  contrasted  with  the  nearly  linear 
dropoff  of  the  forward  current  peak.  The  slight  current  peak  for  the  forward- 
bias  maximum  currents  results  from  a  large  voltage  overshoot  in  the  turn-on 
transient.  The  unfilled  triangle  points  of  the  reverse-bias  curves  were 
calculated  with  initial  space-charge-free  conditions;  the  filled-in  circle 
points  were  calculated  with  stored  charge  resulting  from  a  preceding  forward 
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half  cycle.  As  shown  previously,3  stored  charges  generally  reduce  breakdown 
voltages  and  currents.  The  exceptional  point  at  the  highest  frequency  may  be 
due  to  the  fact  that  the  same  stored  charge  was  assumed  for  the  initial  condi¬ 
tions  for  all  frequencies. 


Figure  7.  Overvoltage  as  a  function 
of  rate  of  voltage  rise  for  data  for 
figure  5.  Assumed  breakdown  voltages 
are  10,  11,  or  12  V.  Slope  of  straight 
line  is  0.58. 


Figure  8.  Comparison  of  forward  and 
reverse  current  peaks  as  a  function 
of  frequency  for  0.  2-jjm  diode.  Pa¬ 
rameter  is  maximum  applied  voltage. 
Empty  circles  and  squares  are  for 
forward  bias,  with  and  without  ava¬ 
lanching,  respectively.  Filled  cir¬ 
cles  and  triangles  indicate  reverse 


respectively. 
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6.  FORWARD  AVALANCHING 

Avalanching  will  also  occur  in  forward-biased  diodes5  with  sufficient 
applied  voltages.  This  is  shown  for  the  100-V  forward-bias  curves  in  fig¬ 
ure  8.  No  difference  is  found  in  dc  calculations  with  or  without  avalanche 
multiplication  since  the  steady-state  voltage  is  too  low.  Comparisons  of 
current  and  voltage  forward  transients  at  160  GHz  are  shown  in  figure  9.  An 
increase  in  current  due  to  avalanching  is  evident  at  about  5  or  6  V;  this  is 
just  about  half  the  reverse-bias  avalanche  breakdown  of  11  V.  This  ratio  is 
typical  of  many  calculations  and  results  from  the  high  forward  (unmultiplied) 
current. 

The  turn-on  transient  of  the  avalanching  forward-current  pulse  of  figure  9 
is  shown  in  figure  10.  The  solid  curves  show  the  doping  profile  of  the  p-i-n 
diode.  Electrons  (wide-dashed  curves)  move  to  the  right,  and  holes  (narrow- 
dashed  curves)  move  to  the  left.  At  the  saturation  velocity,  carriers  move 
100  nm  (0.1  tin)  in  1  ps.  For  reasons  of  clarity,  charge  distributions  during 
the  current  decay  are  not  shown  in  figure  10.  In  this  case  the  current  decay 
is  so  rapid  that  little  charge  redistribution  can  occur.  Since  the  total 
current  must  be  constant  in  space  across  the  diode,  displacement  current  must 
flow  in  the  region  devoid  of  enough  carriers  to  provide  sufficient  conduction 
current.  The  result  is  that  the  field  increases  initially  in  the  center  of 
the  diode.  Field  distributions  during  the  turn-on  transient  of  figure  10  are 
shown  in  figure  11.  Avalanching  occurs  at  fields  about  4  x  105  V/cm,  as  may 
be  seen  from  figures  9  and  1 1 .  This  is  about  one-half  the  reverse-bias 
breakdown  field  for  this  diode. 


& 
v« «  ■ 
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sumption  that  the  thermal  injection  current3  was  the  equilibrium  value  for  the 
indicated  temperature.  For  applied  pulse  lengths  less  than  the  diode  thermal 
time  constant  (of  the  order  of  a  microsecond),  there  is  little  heat  conducted 
and  the  thermal  injection  current  should  correspond  to  a  lower  temperature. 
Figure  13  shows  the  results  for  calculations  made  with  thermal  injection 
currents  corresponding  to  the  100  and  200  K  less  than  the  bulk  temperature. 
Lower  injection  currents  lower  the  cutoff  frequency  but  increase  the  current 
maxima  at  lower  frequencies. 

Further  calculations  of  the  temperature  dependence  of  microwave  current 
peaks  were  made  for  voltage  peaks  of  60  V  applied  to  the  same  diode.  A  second 
difference  in  parameters  was  chosen  to  illustrate  the  effect  of  diode  area 
upon  the  current  peaks.  In  figure  14  the  area  was  chosen  to  be  1  *  10~6  cm2, 
rather  than  the  1  x  10-5  cm2  used  in  figure  12.  The  series  resistance  was 
50  and  the  shunt  capacitance  was  0.2  pF  in  both  cases.  A  comparison  of 
figures  12  and  14  shows  that  the  peak  current-density  maxima  for  the  smaller 
area  diode  are  greater  than  for  the  larger  area  diode,  even  though  the  applied 
voltage  is  less  (see  fig.  5  for  the  variation  of  current  maxima  with  applied 
voltage ) . 


—  Jo  (T-  100) 
--  Jo  (T  -  200) 


FREQUENCY  (GHz) 

Figure  12.  Maximum  current  densities 
as  a  function  of  frequency  at  speci¬ 
fied  temperatures.  Forward-bias  data 
are  shown  as  open  circles  and  reverse- 
bias  data  as  filled  squares. 


FREQUENCY  (GHz) 

Figure  13.  Effect  of  thermal  injec¬ 
tion  currents  (JQ)  upon  maximum  cur¬ 
rent  densities.  Solid  lines  without 
points  are  fully  isothermal  data  of 
figure  12. 
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Figure  14.  Maximum  current  densities 
as  a  function  of  frequency  at  labeled 
temperatures.  Maximum  of  voltage  is 
60  V  and  diode  area  is  1  x  1 0 cm2 . 


Higher  transient  current-density  peaks  for  smaller  area  diodes  result  from 
the  reduced  voltage  drop  in  the  external  resistor.  This  follows  from  the 
lower  total  diode  current  (current  density  times  area)  and  consequently  a 
higher  voltage  across  the  diode.  The  lower  cutoff  frequency  shown  in  figure 
14  compared  to  figure  12  is  a  result  of  the  lower  applied  voltage  (compare 
fig.  5).  There  is  little  difference  in  the  300  and  400  K  curves  in  figure  14; 
this  is  because  the  higher  breakdown  voltage  at  400  K  results  in  higher  cur¬ 
rent  peaks  (compare  fig.  6),  partially  offsetting  the  reduced  ionization 
coefficient  and  saturation  velocity. 

Although  the  initial  temperature  was  assumed  constant  across  the  diode, 
calculations  were  also  made  with  the  option  for  the  temperature  to  increase 
due  to  power  being  dissipated  in  the  diode.  Figure  15  shows  the  diode  voltage 
and  current  density  as  a  function  of  time  for  one  point  (4.8  GHz)  on  the  300  1C 
curve  of  figure  14.  The  breakdown  voltage  of  this  diode  is  extrapolated  to  be 
40  V.  The  maximum  dissipated  power  during  the  pulse  is  just  under  80  W,  and 
the  average  over  the  reverse  half  cycle  is  about  5  W.  Figure  16  shows  the 
average  diode  temperature  as  a  function  of  time  for  this  simulation.  Although 
the  increase  in  average  temperature  is  just  over  3  K  for  this  cycle,  the  maxi¬ 
mum  rate  of  temperature  rise  at  the  current  peak  is  450  K  per  nanosecond.  The 
average  rate  over  the  entire  cycle  is  about  15  K  per  nanosecond.  The  distri¬ 
bution  of  temperature  across  the  diode  at  selected  times  is  shown  in  figure 
17,  which  shows  that  the  temperature  peaks  at  or  near  each  end  of  the  diode. 
These  peaks  are  a  result  of  the  high  space-charge-induced  electric  fields  in 
these  regions.  In  reality,  the  fields  and  temperatures  should  drop  sharply 
from  the  peak  at  each  boundary.  Their  failure  to  do  so  resulted  from  too-low 
doping  levels  assumed  for  these  regions.  These  lower  doping  levels  were  used 
to  allow  calculations  to  be  made  at  lower  currents  and  at  lower  costs. 
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Figure  15.  Reverse  current  density 
and  voltage  as  a  function  of  time 
for  a  1-pm  diode  at  4.8  GHz. 


Figure  16.  Average  diode  temperature 
and  rate  of  temperature  rise  as  a 
function  of  time. 
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Figure  17.  Distribution  of  tempera¬ 
ture  across  diode  at  specified  times. 

The  temperature  distribution  during  the  forward  half  cycle  has  a  peak  in 
the  center  of  the  intrinsic  region.  The  temperature  distribution  is  almost 
identical  in  shape  to  the  field  distribution  shown  in  figure  11.  Calculations 
have  been  made  for  forward  half  cycles  with  an  extrapolated  temperature  dis¬ 
tribution  (essentially  parabolic)  and  for  the  same  constant  temperature  as  the 
average  of  the  dist'ibuted  curve.  The  results  of  the  calculations  differed 
minimally.  Previously,3  similar  results  were  obtained  for  reverse-bias  calcu¬ 
lations.  Since  the  center  of  the  diode  heats  faster  than  the  ends  during  the 
forward  conduction,  and  the  ends  heat  faster  during  reverse  conduction,  the 
final  temperature  distribution  should  be  fairly  uniform  if  the  forward  power 
and  reverse  power  are  roughly  equal. 
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The  rate  of  temperature  rise  for  a  given  diode  at  a  given  applied  voltage 
was  found  to  be  nearly  independent  of  temperature  for  these  calculations.  The 
independence  results  from  the  opposing  effects  of  higher  breakdown  voltage  but 
lower  current  peaks  as  the  temperature  increases.  Knowing  the  increase  in 
temperature  for  both  the  forward  and  reverse  single  half  cycles,  one  can 
readily  approximate  the  total  number  of  cycles  to  obtain  any  given  tempera¬ 
ture.  Thus,  the  time  to  reach  that  temperature  is  found  as  a  function  of  the 
impressed  voltage  amplitude. 

Previous  calculations3  of  single  reverse-bias  square  pulses  indicate  that 
thermal  second  breakdown  occurs  at  roughly  600  K.  It  is  postulated  that  a 
negative  resistance  leads  to  current  constriction  and  burnout.  For  an  applied 
microwave  pulse  at  the  higher  frequencies,  the  time  required  for  current 
contraction  may  have  a  large  effect.  Two-dimensional  calculations  will  be 
required  to  study  this  process. 

8.  CONCLUSION 

Both  forward  and  reverse  characteristics  have  been  calculated  for  micro¬ 
wave  (half -cycle)  voltage  waveforms.  The  cutoff  frequency  for  forward  current 
flow  has  been  shown  to  depend  on  diode  width.  The  reverse  current  flow  has 
been  shown  to  be  important  above  the  diode  breakdown  voltage  and  to  be  sharply 
peaked  in  frequency.  This  shows  that  microwave  burnout  of  diodes  should  be 
strongly  frequency  dependent  at  higher  frequencies  and  short  pulse  lengths. 
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